Measurement of the flux and zenith-angle distribution of upward through-going 

muons in Kamiokande II+III 
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\Q . The flux of upward through-going muons of minimum (mean) threshold energy >1.6 (3.0) GeV 

is measured, based on a total of 372 events observed by the Kamiokande II+III detector dur- 
■ ing 2456 detector live days. The observed muon flux was $ (, s = (1.94 ± 0.10(stat.)^Q Q 6 (sys.)) x 

10 _13 cm _2 s -1 sr -1 , which is compared to an expected value of $>theo = (2.46 ± 0.54(theo.)) x 
' 10 _13 cm _2 s _1 sr _1 . The observation is in agreement with the prediction within the errors. The 

^ i zenith angle dependence of the observed upward through-going muons supports the previous indica- 

tion of neutrino oscillations made by Kamiokande using sub- and multi-GeV atmospheric neutrino 
i-^ ' events. 

> : 

In 1988, the Kamiokande group observed a smaller (u u + Vy)j(v e + v e ) ratio than the expected value in sub-GeV 



fully-contained atmospheric neutrino events [[L||2|] . This discrepancy was later confirmed by other experiments IMB-3 
HQ and Soudan II ||, while some other experiments (Frejus || and Nusex 0) did not observe such discrepancy 
within the experimental errors. 

Recently, a measurement of this ratio using multi-GeV atmospheric neutrino events was made by the Kamiokande 
experiment || . The observed (y^ +i7 /1 )/(i/ e +V e ) ratio in this energy range was also significantly smaller than expected 
and depended strongly on the zenith angle. 

Having observed the discrepancy in both sub- and multi-GeV atmospheric neutrino events, the next logical step 
|p|, ^0|JTl]Jl^| is to investigate events that are produced by neutrinos of even higher energies. Energetic atmospheric 
or traveling through Earth interacts with rock layers surrounding the detector and produce muons via charged- 
current interactions. Though neutrino- induced downward-going muons hit the detector, these downward-going muons 
are difficult to differentiate from cosmic-ray muons. On the contrary, upward-going muons are considered neutrino- 
induced, since upward-going cosmic-ray muons range out in Earth. Especially, those energetic enough to pass through 
the detector are defined as "upward through-going muons" . The mean energy of the parent neutrino producing them 
is approximately 100 GeV. 

The experimental site is located 1000 m underground in the Kamioka mine, Gifu prefecture, Japan. The 
Kamiokande-II detector is a cylindrical water Cherenkov calorimeter. The inner detector is 14.4 m in diameter 
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x 13.1 m in height, which contains 2142 metric tons of purified water. A total of 948 20-inch-in-diameter photomul- 
tiplicr tubes (PMTs) are mounted on the tank walls and covers 20% of the inner tank surface. The detector also has 
an optically isolated 4ir solid-angle anticountcr which is also a 1500-metric-ton water Cherenkov detector with 123 
20-inch-in-diameter PMTs and is divided into top, barrel and bottom parts. The anticounter tags incoming and/or 
outgoing muons and to shield 7-rays and neutrons from the surrounding rock. Charge and relative timing information 
from each hit PMT is recorded for the event reconstruction. The triggering efficiency for a muon having momenta 
more than 200 MeV/c is ~100%. The nominal detector effective area for upward-going muons is roughly 150 m 2 . 
The deviation between the reconstructed track direction and the real muon direction (A9 rec ) is studied by using 
Monte Carlo events and is estimated to be 2.1° [||. See 0,0] for more details of the Kamiokande-II detector. The 
new phase of the experiment started in December 1990, after replacement of the entire electronics, replacement of the 
dead PMTs and installation of an aluminized mylar cone around each PMT. The description of the Kamiokande-III 
detector is found elsewhere [jlj|. 

The data sets used in this analysis are taken from December 1985 to April 1990 in Kamiokande-II and from 
December 1990 to May 1995 in Kamiokande-III corresponding to 2456-day detector livetime. A total of 2.2 x 10 8 
events were recorded during these data-taking periods. 

In Kamiokande-III, when the following three requirements are satisfied, it is considered that a muon hits one of the 
three anticounter parts: (1) total number of photoelectrons (p.e.) in a part of the anticounter > 20 p.e. (2) number of 
hit PMT in the same part of the anticounter > 5 (3) Maximum number of photoelectrons/PMT in an event in the same 
part of the anticounter > 5 p.e. By requiring the two parts or the barrel part hit in the anticounter, through-going 
muons can be selected. A selection criterion that the total number of photoelectrons of an event should be larger than 
6000 p.e. is imposed to ensure the minimum muon track length of 7 m (corresponding to minimum (mean) muon 
threshold energy of 1.6 (3.0) GeV). An upper limit is also set at 40000 p.e. (corresponding to ~ 30 m track length 
for a relativistic muon in the detector) to remove multiple muons or muons accompanied with bremsstrahlung, pair 
production and/or hadronic showers. The event reduction procedures in Kamiokande-II is described elsewhere ||. 

To eliminate abundant downward-going cosmic-ray muons (0.37 Hz), the events satisfying cos0 > —0.04 are cut. 
where is the zenith angle with cos = 1 corresponding to downward-going events. The detection efficiency for 
upward through-going muons is estimated by a Monte Carlo simulation. Using the upward/downward symmetry of 
the detector configuration, the validity of the Monte Carlo program is checked by real cosmic-ray downward through- 
going muons. The average detection efficiency is estimated to be ~ 97%. The details of the selection criteria are 
described in [l3| . Finally, 184 and 188 upward through-going muon events are observed during Kamiokande-II PJT^] 
(1124 days) and Kamiokande-III jy| (1332 days) experimental periods, respectively. 

In this analysis, the combination of the Bartol atmospheric neutrino flux model |l5f| , the parton distribution functions 
of GRV94DIS [^6| and the Lohmann's muon energy loss formula in the standard rock jL7| are employed, when one 
analytically calculates the standard expected muon flux. Also a Monte Carlo technique is introduced in order to 
estimate the directionality between a muon and the parent neutrino and it is found that J AO^ + AO^ui is 4.1° for 
the atmospheric neutrino energy spectrum, where A9 V ^ stands for /1 production angle relative to v direction in the 
neutrino-nucleon interaction and A9 mu i for deflection by multiple Coulomb scatterings in the rock. To estimate the 
model-dependent uncertainties (±10 % for the absolute flux normalization and —3.6% to +1.5 % for the bin-by-bin 
shape difference in zenith-angle distribution) of the expected muon flux, the combination of another atmospheric 
neutrino flux model calculated by Honda et al. (ll| and another parton distribution functions of CTEQ3M are 
also considered. 

As a result, the expected muon flux $ t / leo is calculated to be (2.46±0.54(theo.)) x 10~ 13 cm -2 s _1 sr _1 (cos0<— 0.04), 
where the estimated theoretical uncertainties are described in Table || The dominant one originates from the absolute 
normalization uncertainty in the neutrino flux which is estimated to be approximately ±20 % [p^| , ^8|j2C| l above several 
GeV. 

Given the detector live time, T, the effective area for upward through-going muons, S(@), and the detection 
efficiency, e(0), the upward through-going muon flux is calculated by the formula: 
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£2(e j )T2 + e 3 (e j )T 3 J 5(9^)2(1-0.04)77 



where suffix 2 (3) represents Kamiokande-II (-III), the suffix j stands for each event number, 2(1 — 0.04)7r is the total 
solid angle covered by the detector for upward through-going muons, N corresponds to the total number of observed 
muon events (372 events). The angular dependence of the detection efficiency for each experimental period is found 
in H and [fl3"[ , respectively. Conceivable experimental systematic errors are summarized in Table |[ 

Taking these experimental systematic errors into account, the observed upward through-going muon flux is: 



$ obs = (1.94±0.10(stat.)To;o 6 ( s y s -)) x 10~ ia cnr W" 1 . 

At present, three experiments other than Kamiokande reported the measurement of the upward through-going 
muon flux and their results are summarized in Table |fi|, together with this measurement. The observed flux by each 
experiment is in agreement with the expectation within the errors. 

The zenith angle distribution of the observed flux {d<& / dVtY obs is shown in Fig. [jj The shape of the distribution is 
not well represented by the theory (% 2 /degrees of freedom = 21.3/9.) 

A neutrino oscillation hypothesis is then tested using the zenith angle distribution. The expected flux for a given 
set of Am 2 and sin 2 29 is calculated and the same binning is applied as data. 

To test the validity of the oscillation hypothesis, a x 2 is calculated as 

\ 2 = min 

where <J s tat,i {csys,i) is the statistical (experimental systematic) error in (d<fr/dQy obs in the i-th bin, a is an absolute 
normalization factor of the expected flux. Based on uncorrelated systematic errors in Table [| added in quadrature, 
we estimate a sys ,i to be ±18.5 % for — 0.1<cos8<— 0.04 and ±1.2 to ±3.8 % for —1 <cos0 < —0.1, respectively. 
Although the nature of the uncertainty in a is unknown, it is presumed that a obeys a Gaussian distribution with 
one standard deviation error a a . The absolute flux normalization error a a is estimated to be ±22 % by adding in 
quadrature the correlated experimental errors and theoretical uncertainties in Table || Then, the minimium X 2 (Xmm) 
on the Am 2 — sin 2 29 plane is searched for. 

We note here that the allowed region given by Kamiokande sub- and multi-GeV contained event analysis j^] for 
v e <-* oscillations has been already excluded by the CHOOZ experiment [p3[ , suggesting <-» v T oscillations for 
the contained and upward-going muon events. Assuming «-» v T oscillations, Xmin( = 12.2) occurs in the unphysical 
region at (sin 2 20, Am 2 ) = (1.35, 2.0xlO _3 eV 2 ) anda=1.02. On the other hand, if we bound the oscillation parameters 
in the physical region, xL™(= 12 - 8 ) takes P lace at (sin 2 26», Am 2 ) = (1.00,3.2 x 10~ 3 cV 2 ) and a=1.00. For the null 
oscillation case, we obtain Xmin °f 21.3 at a=0.77 (probability of statistical fluctuation: 1%). The zenith angle 
distribution of a{d<& / d£l) l osc for the best fit parameters in the physical region is shown in Fig. |l| together with the 
data. Figure |^ shows the allowed region contours on the (sin 2 29, Am 2 ) plane for <-» v T oscillations. As Xmin f° r 
<-> v T oscillations falls down in the unphysical region, the contours are drawn according to the prescription for 
bounded physical regions given in Ref. If we replace the Bartol neutrino flux [|l5| by the Honda's jl8) and/or the 
GRV94DIS parton distribution functions j[6| by the CTEQ3M's the allowed region contours are similar to those 
presented in Fig. |2| Consequently, we find that the zenith angle dependence is in favor of the <-> v T oscillation 
hypothesis and supports the Kamiokande sub- and multi-GeV contained event analysis Q. 

Combining this work with the previous analysis || of the Kamiokande sub- and multi-GeV contained events, we 
draw an allowed region contour on the (sin 2 29, Am 2 ) plane, as is shown in Fig. 0. The best fit point (sin 2 29, Am 2 ) — 
(0.95, 1.3 x 10- 2 eV 2 ) is obtained for Vp, *-> v T oscillations. 

In conclusion, based on 372 upward through-going muon events during 2456 detector live days, the flux of the upward 
through-going muons (>1.6 GeV) produced by atmospheric muon neutrinos in the rock layers surrounding the detector 
is measured with the Kamiokande II+III detector: $ ohs = (1.94 ± 0.10(stat.)^p'pg(sys.)) x 10 _13 cm _2 s _1 sr _1 . This 
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is compared with the expected flux calculation of $theo = (2.46 ± 0.54(theo.)) x 10 _13 cm _2 s _1 sr _1 . The observed 
upward through-going muon flux is in agreement with the expected flux within the relatively large uncertainties in 
the theoretical calculations. We find that the zenith angle dependence of the upward through-going muons does not 
agree with the theoretical expectation (probability of statistical fluctuation: 1%), yet is favored by the <-> v T 
oscillation assumption. This result supports the indication of neutrino oscillations given by the analysis of the sub- 
and multi-GeV atmospheric neutrino events by Kamiokande. 
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FIG. 1. Zenith angle distribution of upward through-going muon flux observed in Kamiokande II+III. Inner (outer) error bars 
indicate statistical (uncorrelated experimental systematic + statistical added in quadrature) errors. The solid histograms show 
the expected (Bartol + GRV94DIS + Lohmann's ionization formula) upward through-going muon flux for the null neutrino oscil- 
lation case. Also shown is the expected flux (dashed) assuming the best fit parameters at (sin 2 28, Am 2 ) = (1.00, 3.2 x lCP 3 eV 2 ), 
a=1.00 in the physical region for the <-> v T oscillation case . 
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FIG. 2. The allowed regions at 90 % C.L., obtained from the Kamiokande II+III upward through-going muon analysis(up 
fi) and from the combined analysis (Kamiokande contained M + up n), drawn on the (sin 2 29, Am 2 ) plane for <-> v T 
oscillations. Also shown are the excluded region by CDHSW |25| and the allowed region by the Kamiokande contained event 
analysis (jl. The right-hand side of the contours are allowed (excluded) regions by Kamiokande (CDHSW). The asterisk 
(cross) indicates the best fit point for up n (for Kamiokande contained + up /i) at (sin 2 29, Am 2 ) — (1.00,3.2 x 10 _3 eV 2 ) (at 
(sin 2 29, Am 2 ) = (0.95, 1.3 x 10~ 2 eV 2 )). The 95 % C.L. contour (up ^) is presented as well. 
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Error source 


Error (%) 


uncertainty in A8 rec 


±l a 


cosO = —0.04 cut dependence 


+ 18.2, -4.4 a 


cosmic ray /i contamination 


-10 a 


detection efficiency 


+2.8, -8.3 a 




+0.6, -1.2 b 


7 m track length cut 


±0.5 C 


live time 


±0.4 C 


effective area 


±0.3 C 


PMT gain 


<l c 


water transparency 


<l c 


chemical component in the rock 


«l d 


v flux normalization?^ 


±20 d 


theoretical model dependence 


±10 d 




-3.6 to +1.5 C 



TABLE I. List of experimental systematic errors and theoretical uncertainties in the flux measurement, a: experimental 
uncorrelated systematic error specific in the most horizontal bin — 0.1<cosO<— 0.04, b: experimental uncorrelated systematic 
error specific in the bins —1 <cos0 < —0.1, c: experimental correlated systematic error, d: theoretical correlated uncertainty, 
e: theoretical uncorrelated uncertainty. 
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Experiment 


Minimum E^(GeV) 


Observed 


Expected 


Kamiokande 


1.6 


It n A I n -in/ j_ j_ \-t-0.07/ w 

(1.94 ± 0.10(stat.)I g 06 (sys.)) 
x 10~ 13 cm _2 s sr -1 


(2.46 ± 0.54(theo.)) 

x 10 _13 cm~ 2 s _1 sr _1 


1MB p) 


~1.4 


532±23(stat.) events 


516±103(theo.) events 


Baksan 


1.0 


559 events 


580 events 


MACRO |J 


1.0 


277 ± 17(stat.) ± 22(sys.)events 


371 ± 63(theo.)events 


TABLE II. 


Summary of upward-going muon flux 


measurement. E M represents the muon energy. 


Kamiokande, Baksan and 



Macro employ the Bartol neutrino flux Jl5j, while 1MB adopts the Lee and Koh neutrino flux j2lj (the Volkova flux |2j) for 
neutrino energies lower (higher) than 15 GeV. 
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